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1492keshav @ gmail.com inadequate thermal performance and require mechanical devices to bring thermal 
comfort. Hence, it is necessary to know the thermal performance of the RCC roof 
structure, which is quite commonly available. This paper concentrates on developing a 
Thiagarajar College of numerical model of a Reinforced Concrete (RCC) roof for analyzing its thermal per- 


formance. Numerical Simulations have been performed for unsteady conditions with 
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I. INTRODUCTION 


Climate change, caused by the release of greenhouse 
gases (mainly carbon dioxide) into the atmosphere, has 
been recognized as one of the greatest threats of the 
21st century. The total amount of energy used by 
commercial buildings has risen significantly since the 
1980s, reflecting a 50% growth in the total amount of 
office space available and a 33% increase in energy 


time-varying solar heat flux along with grid independence tests for accurate analysis. 


Index Terms— RCC Roof, Thermal Performance, Numerical Simulation 


50 percent of lime mortar by volume is laid over the 
RC slab. The concrete is well rammed so that the 
thickness reduced, by wooden hand beaters. The sur- 
face 1s cured by sprinkling lime water. When the brick 
bat is set a course of flat tile 1s laid in lime mortar (1:1 
Y). This type of roof is conventionally used in most of 
the residences in hot humid climates in Tamil Nadu, 
India hence it 1s taken as a reference roof for study. 
The roof section is schematically represented in Fig. 1. 


The thermo physical properties [5] of roof materials are 
shown in Table 1. 


consumption per square foot of space [1]. With the 
convergence of urbanization, globalization and a rapid- 
ly changing and expanding economy, India is expe- 
riencing a rapid spurt in building construction across a 
range of city activities and socio-economic spectrum, 
increasing consumption of building materials such as 
glass, cement, metals and ceramics. Maximum con- 
sumption of these energy materials is a reason for envi- 
ronmental degradation. Buildings are the dominant 
energy consumers in modern cities accounting up to 
40% energy consumption. There is over 50% saving 
potential in the building sector and thus it is considered 
as a potential sector to meet the challenges of global 
energy and climate change [1]. Roofs are a critical part 
of the building envelopes that are highly susceptible to 
solar radiation and other environmental changes, there- 
by, influencing the indoor comfort conditions for the 
occupants [2]. In developing countries, masonry houses 
with RCC roofs are popular owing to their pest resis- 
tance, natural calamities resistance, availability and 
cost effectiveness of concrete ingredients [3]. Depend- 
ing on the weather, the thermal load enters into a build- 
ing in three major ways such as penetration of direct 
beam sunlight, conduction of heat through roofs, infil- 
tration of outside air [4]. 





Fig. 1 RCC roof section 


Table 1 Properties of materials 


MATERIAL | DENSITY SPECIFIC THERMAL 
; | HEAT CONDUCTIVITY 
(kg /m”) 
(J/kg-K) (W/mK) 
CLAY 1845 1089 1.036 
CON- 2300 1130 1.27 
CRETE 
CEMENT 1900 1130 0.3024 


Ill. GEOMETRIC MODELLING 


The roof model is made in GAMBIT software. The 
model is made in the form of 2-dimensional structure 
of 180x100 mm. The leftmost portion of model is de- 

fined as “outside” of the roof. The rightmost portion of 


Il. DESCRIPTION OF RCC ROOF 





Flat clay tile roof 1s selected for modeling and simu- 
lation. Horizontal slabs of steel reinforced concrete, 
typically between 100 and 500 millimeters thick, are 
most often used to construct floors and ceilings. A 10 
cm thick layer of brick-bat concrete is laid, consisting 
of 3 parts of brick-bats, | part of gravel and sand, and 
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the model is defined as “inside” of the roof. All the 

four portions of the roof, of different materials, are 

defined as separate faces. Fig. 2 shows the meshed 
model. 





Fig. 2 Meshed Model of roof section 


IV. GRID INDEPENDENCE TEST 


The meshing was made with different interval sizes, 
ranging from Q.1 to 0.9. Each and every model was 
subjected to steady state conditions for the analysis 
using constant heat flux of 700W/m” defined on out- 
side of the roof and convective heat transfer at inside 
portion of the roof. Fig. 3&4 show the temperature 
plots. From the Fig. 3&4, it can be concluded that, the 
curve seems to saturate around interval size of 0.8.The 
temperature difference matches closely for interval size 
of 0.7 and 0.8.The time taken for iteration is lesser for 
interval size 0.8 compared to 0.7. Thus 0.8 1s chosen as 
optimum grid size. 
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Fig. 3 Temperature Vs Interval Size plot 
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Fig. 4 Temperature Difference Vs Interval Size plot 
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IV. NUMERICAL SIMULATION OF RCC ROOF 


In a normal day, the solar heat flux is never a con- 
stant quantity. It is estimated that average heat flux on 
a normal sunny day, increases in magnitude from 
morning to noon, reaches the peak value at afternoon 
and then decreases during evenings. It 1s observed that 
in a normal day, influence of solar heat prevails for 8 
hours [6]. The residents feel comfort during the period 
of the measurement, in which the Mean Room Tem- 
perature ranged 27~32°C, the room air temperature 
28~32°C [7]. Thus, the room air temperature is taken 
as 303K. Unsteady conditions have been taken into 
account. The convection effect is considered around 
the roof edges with heat transfer coefficient assumed as 
10W/m? K. It is necessary to use time dependent input 
as boundary condition. Time varying heat flux input is 
obtained with the help of User Defined Function 
(UDF). Using [8] the UDF was developed according to 
the required condition. The iteration is performed for a 
period of 8 hours. The initial temperature of the roof is 
assumed as 300K. The convergence criterion 1s set as 
le-08. Fig. 5 shows the applied input heat flux, as a 
plot. 
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Fig. 5 Input Heat Flux plot 


The obtained temperature contour is shown in Fig. 
6. It is observed that, maximum temperature is 332K 
and minimum temperature is 306K after 8 hours. This 
drop occurs due to varying conductivities of the differ- 
ent materials in the roof section. During simulation, 
the average outside and inside roof temperature data 
are obtained and plotted in a graph as in Fig. 8. 
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Fig. 6 Temperature Contours 
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Fig. 7 Enthalpy Contours 
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Fig. 8 Temperature Vs Hours plot 


VI. VALIDATION 


The actual conditions are determined using non- 
contact temperature measuring equipment (Infrared 
camera). The thermal image is captured in the after- 
noon, when the solar radiation is intense. The captured 
image is displayed above in Fig. 9. The outside tem- 
perature obtained from the simulation, is in agreement 
with the thermal image result and literature [9]. 


65.9°C 





Fig. 9 Thermal Image of outside portion of the roof 
V. RESULT AND DISSCUSSION 


The temperature contours (Fig. 6) show the tempera- 
ture distribution along the roof. It is observed that tem- 
perature drop of 31.6 K occurs along the roof at peak 
condition, when solar intensity is maximum. This 
means the input energy is stored in the roof. The en- 
thalpy variations show that, enthalpy, as a function of 
internal energy, is dependent on the specific heat ca- 
pacity of the material. Another significant observation 
made from the plot is that, flat clay tile has higher en- 
thalpy compared to other materials. The cement layer 
has lower conductivity. It does not transfer energy 
from concrete quickly. From the Fig. 8, the inside 
temperature increases linearly with hours. The outside 
temperature varies in a non-uniform manner as it de- 
pends on solar heat flux .Further, it is found to increase 
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till 6 hours and then decrease in accordance to heat 
flux. 


VI. CONCLUSION 


The thermal performance of the RCC roof section is 
numerically analyzed. The grid independence test was 
conducted and the optimum grid size has been deter- 
mined. The temperature distribution across the roof has 
been obtained. This paper will prove useful for estima- 
tion of heat load of RCC structure. Currently, a lot of 
advancement techniques are springing up for the reduc- 
tion of heat loads by varying the passive design. This 
paper can provide basic information for such ideas. 
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